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Preface 


This  computer  study  Is  conducted  to  determine  the  effects  of 
concentrated  mass  combinations,  simulating  external  stores,  on  the 
flutter  characteristics  of  a non-uniform  wing.  This  approach,  I 
believe,  provides  a rapid,  versatile  method  of  predicting  flutter 
problems  with  aeroelastlc  structures  during  the  preliminary  design 
process,  where  flutter  problems  can  most  easily  and  economically  be 
remedied.  This  study  should  be  of  Interest  to  the  aerodynamlclst  and 
the  aeroelastlclan,  as  flutter  Is  a concern  of  both. 

An  advanced  version  (Level  16.0)  of  the  NASA  Structural  Analysis 
(NASTRAN)  computer  system  Is  used  to  perform  the  flutter  analysis.  A 
finite  element  model  of  a strength  optimized  wing  Is  used  In  conjunc- 
tion with  doublet- lattice  aerodynamics. 

There  are  two  limitations  of  this  study.  One  Is  that  the 
aerodynamic  Interaction  between  the  structure  and  attached  stores 
has  not  been  Included  In  the  computer  model.  The  second  Is  the 
availability  of  only  doublet-lattice  aerodynamics  In  NASTRAN,  which 
restricts  this  method  to  subsonic  lifting  surfaces  undergoing  small 
sinusoidal  motions. 

I would  like  to  thank  my  thesis  advisor.  Major  Franklin  Eastep, 
and  my  faculty  committee  for  their  guidance  and  assistance.  In  addition, 
I would  like  to  thank  Captain  Gene  Hemmlg  for  his  assistance  In  the 
application  of  NASTRAN.  Finally,  I would  like  to  thank  my  wife  Cheryl 
for  her  pat  ^nd  support  during  this  project. 


Van  C.  Sherrer 
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Abstract 

^j| 

' A computer  study  of  the  effects  of  external  stores  simulated  by 
limped  masses  was  conducted  with  a finite  element,  cantilever,  non- 
uniform  wing  model.  The  NASTRAN  (Level  16.0)  computer  program  flutter 
format  was  used  to  obtain  flutter  speeds  and  frequencies.  Mass  balancing 
with  a single  concentrated  mass  caused  a reduction  In  flutter  speed  as 
the  mass  was  moved  chordwlse  toward  the  trailing  edge  and  spanwlse 
toward  the  wing  tip.  Flutter  speeds  and  frequencies  of  a 100  lb  and  a 
200  lb  store,  simulated  by  two  equal  masses,  were  compared  to  an  equiv- 
alent concentrated  mass  at  the  store  center  of  gravity.  The  stores 
consistently  raised  the  flutter  frequency  over  that  of  the  single  mass, 
but  flutter  speed  results  were  not  conclusive. 
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THE  EFFECTS  OF  EXTERNAL  STORES  ON  THE 


FLUTTER  OF  A NON-UNIFORM  CANTILEVER  WING 

I.  Introduction 

The  purpose  of  this  computer  study  Is  to  Investigate  the  effects 
of  external  stores,  simulated  by  concentrated  mass  combinations,  on  the 
flutter  characteristics  of  a non-uniform  cantilever  wing.  A non-uniform 
wing  is  defined  In  this  study  as  a wing  with  varying  stiffness  and  mass 
properties  along  the  span  of  the  wing.  This  type  of  analysis  will  be 
useful  In  preliminary  design  studies  as  a method  of  rapidly  estimating 
the  flutter  speed  of  an  aeroelastlc  structure  as  It  evolves  during  the 
design  process.  Flutter  speed  determination  Is  required  when  the  mass 
and/or  stiffness  of  the  structure  Is  Increased  or  decreased  as  a result 
of  design  changes,  and  during  stores  flutter  clearance  studies.  If  the 
predicted  flutter  speeds  are  lower  than  desired,  then  changes  can  be 
Incorporated  at  an  early  stage  of  the  design  process  to  raise  the  flutter 
speed  or  eliminate  flutter  altogether.  The  earlier  a design  problem  Is 
detected,  the  less  expensive  It  Is  to  remedy. 

The  flutter  of  wings  with  external  stores  Is  a subject  of  current 
Interest  In  the  Air  Force.  Each  aircraft  that  carries  external  stores 
must  be  cleared  for  Its  particular  flight  envelope  for  each  unique  store 
combination,  or  else  operating  restrictions  must  be  Imposed  on  the  air- 
craft In  order  to  guarantee  a margin  of  safety  from  flutter.  The  number 
of  external  ?tore  combinations  can  reach  astronomical  proportions  with 
certain  aircraft,  thus  the  need  for  a rapid  flutter  prediction  method 
exists. 
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Flutter 


Definitions.  Flutter  Is  described  by  Blspllnghoff  as  the  "dynamic 
Instability  of  an  elastic  body  In  an  alrstream"  (Ref.  2:527).  The  flutter 
speed  Is  defined  as  the  lowest  airspeed,  U^,  at  which  a given  structure 
flying  at  a given  atmospheric  density  and  Mach  number  will  exhibit  sus- 
tained, simple  harmonic  motion.  The  corresponding  circular  frequency, 
u^,  at  which  the  flutter  speed  occurs  Is  the  flutter  frequency.  The 
flutter  speed  represents  a boundary  condition  below  which  stable  oscilla- 
tions of  the  structure  occur.  At  speeds  above  the  flutter  speed,  divergent 
oscillations  with  resultant  structural  failure  can  occur. 

A static  Instability  problem  of  lifting  surfaces,  known  as  divergence, 
should  also  be  mentioned.  The  most  common  case  of  divergence  Is  torsional 
divergence,  which  occurs  when  the  Increment  In  aerodynamic  torsional 
moment  caused  by  an  Increase  In  twist  angle  exceeds  the  elastic  restoring 
torque  of  the  wing  structure.  At  that  point,  the  wing  structure  suffers 
a catastrophic  failure.  Divergence  speed,  however.  Is  usually  higher 
than  flutter  speed  for  swept  back  wings. 

Background.  Historically,  flutter  problems  began  to  appear  as 
airplanes  became  capable  of  higher  and  higher  airspeeds.  To  reach  these 
higher  speeds,  designers  began  to  abandon  biplane  construction  with  Its 
relatively  high  torsional  stiffness  and  aerodynamic  drag  in  favor  of 
lower  drag  monoplane  designs,  which  had  Insufficient  torsional  stiffness 
to  prevent  flutter.  During  the  development  of  monoplane  aircraft,  serious 
research  Into  the  aeroelastlc  problem  began. 

The  basic  two-dimensional  theory  of  flutter  was  originally  devised 
by  Theodore  Theodorsen  (Ref.  12)  in  1934.  He  recognized  that  classical 
flutter  Is  characterized  by  the  Interaction  of  aerodynamic,  elastic,  and 
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inertial  forces.  Additionally,  this  Interaction  Is  usually  the  result  of 
the  coupling  of  two  or  more  degrees  of  freedom.  In  order  to  prevent  this 
coupling,  either  the  stiffness  of  the  wing  Is  Increased  or  the  mass  Is 


redistributed,  or  both. 

In  1940,  Theodorsen  and  Garrick  (Ref.  13)  Investigated  the  effects 
of  various  parameters  on  flutter,  such  as  center  of  gravity  location  of 
a wing  section,  structural  friction,  radius  of  gyration,  and  flutter 
frequency.  They  concluded  that  normally  the  most  important  parameter  Is 
the  center  of  gravity  location.  In  general,  the  further  aft  the  location 
of  the  chordwlse  center  of  gravity,  the  lower  the  flutter  speed. 

In  1945,  Goland  (Ref.  5)  was  able  to  solve  the  bending-torsion 
flutter  problem  exactly  for  a uniform  cantilever  wing  with  constant 
chord  and  uniformly  distributed  mass  and  elastic  properties.  This  solution 
was  then  used  as  a standard  of  comparison  to  estimate  the  accuracy  of  more 
approximate  methods. 

A wind  tunnel  study  of  the  effects  of  concentrated  weights  on  the 
flutter  characteristics  of  a straight  cantilever  wing  model  was  conducted 
by  Runyan  and  Sewall  (Ref.  11)  in  1948.  Their  general  result  was  that  as 
mass  was  moved  chordwlse  from  the  leading  edge  to  the  trailing  edge, 
flutter  speed  was  reduced.  As  the  weight  was  moved  spanwlse  from  root 
to  tip,  there  was  a general  reduction  of  flutter  speed  and  then  an  Increase 
at  the  tip  as  compared  to  the  unweighted  wing.  The  flutter  mode  was 
generally  first  bending  for  inboard  positions  of  the  weights  and  generally 
second  bending  for  weights  at  the  tip. 

Since  the  advent  of  high  speed  computer  technology  and  finite  element 
techniques  for  both  structures  and  aerodynamics,  it  Is  now  possible  to 
solve  the  flutter  problem  rapidly  and  accurately.  Several  computer  programs 
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have  been  developed  specifically  for  the  solution  of  the  flutter  problem 
In  recent  years.  An  existing  structural  analysis  program  has  now  Incor- 
porated the  flutter  solution  In  Its  latest  version.  This  program,  which 
la  used  In  this  study.  Is  described  In  the  next  section. 

NASTRAN 

The  NASA  Structural  Analysis  (NASTRAN)  computer  system  (Level  16.0), 
developed  under  contract  for  NASA  by  a team  composed  of  Computer  Sciences 
Corporation,  the  MacNeal-Schwendler  Corporation,  and  the  Martin  Company, 
Is  used  In  this  study  to  perform  the  flutter  analysis  (Ref.  8).  NASTRAN 
uses  finite  elements  to  model  the  distributed  physical  properties  of  a 
structure.  The  finite  elements  are  Interconnected  at  grid  points  to 
which  loads  are  applied  and  for  which  displacements  are  calculated.  Mass 
properties  are  either  calculated  Internally  as  properties  of  structural 
elements  or  Input  as  properties  of  grid  points.  Constraints  of  various 
kinds  may  be  applied  to  the  grid  points  to  simulate  boundary  conditions 
or  tie  the  movements  of  grid  points  together. 

Three  rigid  formats  of  NASTRAN  are  used  In  this  study.  The  Static 
Analysis  format  Is  used  to  determine  the  elastic  axis  of  the  wing  model. 
The  Normal  Modes  Analysis  format  Is  used  to  determine  the  natural  mode 
shapes  and  frequencies  of  the  wing  vibrating  in  vacuum.  Finally,  the 
Modal  Flutter  Analysis  format  is  used  to  determine  flutter  speed  and 
frequency  for  each  mass  combination  Investigated. 

With  Level  16.0  of  NASTRAN,  the  rigid  format  for  modal  flutter 
analysis  by  the  k-method  was  added.  NASTRAN  now  has  the  capability  to 
generate  aerodynamic  grid  points.  Interpolate  between  the  structure  and 
aerodynamics,  compute  aerodynamic  matrices,  and  solve  the  flutter 
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eqiiatlons  (Ref.  7:17.1-1).  A general  discussion  of  the  theory  employed  by 
NASTRAN  for  the  flutter  analysis  follows  in  the  next  chapter. 


II.  Theory 

Finite  Elements 

The  matrix  methods  of  structural  analysis  used  by  digital  computers 
are  based  on  the  Idea  of  replacing  an  actual  continuous  structure  by  an 
equivalent  model  made  up  of  a finite  number  of  discrete  structure  elements. 

Each  of  these  elements  has  known  elastic  and  Inertial  properties  that  can 
be  expressed  In  matrix  form.  These  matrices,  when  joined  together  In  a 
manner  consistent  with  a set  of  rules  derived  from  the  theory  of  elasti- 
city, yield  the  static  and  dynamic  properties  of  the  actual  continuous 
structure.  Przemlenleckl' s book  provides  an  extensive  treatment  of  the 
theory  of  matrix  structural  analysis  (Ref.  10). 

Three  types  of  finite  elements  are  used  to  model  the  wing  used  In 
this  analysis:  rods,  shear  panels  and  quadrilateral  plate  elements 
(Ref.  8:1. 3-4,  1.3-5).  The  rod  element  Includes  extenslonal  and  torsional 
properties,  but  It  does  not  resist  bending.  The  shear  panel  Is  a two- 
dimensional  element  that  can  resist  tangential  forces  along  Its  edges,  but 
It  cannot  resist  normal  forces.  The  quadrilateral  plate  element  has  both 
Inplane  and  bending  stiffness  with  a solid  homogeneous  cross  section.  | 

These  elements  are  joined  at  node  points  to  form  the  representation  of 
the  wing  box  structure,  as  seen  In  Fig.  1. 

Doublet  Lattice  Method 

The  basic  technique  employed  to  represent  lifting  surfaces  In  NASTRAN 
Is  the  Doublet-Lattice  Method  (DLM)  as  originally  described  by  Albano  and 
Rodden  (Ref.  1:279-285).  The  DLM  Is  an  extension  of  the  Vortex-Lattice 
Method  to  subsonic  unsteady  flow.  The  lifting  surface  Is  modeled  by 
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horseshoe  vortices  that  represent  the  steady  flow  effects  and  acceleration 
potential  doublets  that  represent  the  oscillatory  effects. 


To  apply  DLM,  the  lifting  surface  Is  divided  Into  arrays  of  trape- 
zoidal boxes  whose  sides  lie  parallel  to  the  freestream  so  that  surface 
edges,  fold  lines,  and  hinge  lines  fall  on  the  box  boundaries.  The 
bound  vortex  and  a distribution  of  acceleration  potential  doublets  are 
placed  on  the  quarter-chord  line  of  each  box,  as  In  Fig.  2.  A control 
point  Is  centered  on  the  three-quarter  chord  point  of  each  box.  Then 
the  effects  of  all  vortices  and  doublets  are  summed  for  each  control 
point  to  obtain  the  total  normalwash  at  a control  point. 


DOWNWASH 
COLLOCATION 
POINT 


BOUND  VORTEX 
AND  LINE  OF 
DOUBLETS  (^1^) 

TRAILING 

VORTICES 


Fig.  2.  Lifting  Surface  Idealization 
by  Doublet-Lattice  Method 
(Ref.  4:1A) 


The  integral  equation  that  relates  the  wash  normal  to  a harmonically 
oscillating  surface  to  the  lifting  pressure  Is 


w(PJ) 


Sir 


(P)  KF(PJ,P,k,m)dS 


(1) 
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where 


( 

w 21c 

w ■ normalwash  angle  - — ■ a + i ( — )h 

U c 

w ■ downvash  velocity 
F * any  point  on  the  planform 
Pj  ■ the  point 

Cp  ■ differential  pressure  coefficient 

KF  ■ kernel  function  relating  normalwash  at  Pj  to  the  unit 
pressure  at  P (Ref.  14) 

cu 

k ■ reduced  frequency  ” 
m ~ Mach  number 

h ■ vertical  displacement  of  the  surface 
a > streamwlse  slope  of  the  deformed  surface 
c ■ reference  chord 
S ■ lifting  surface  area 
U ■ freestream  velocity 
u ■ frequency  of  oscillation 

(Ref.  12:78-79) 

This  integral  is  then  discretized  and  put  in  matrix  form  as  used  in 
NA5TRAN: 

f.j)  - t*j,I  tPj)  (2) 

downvash  (normalwash)  at  3/4  chord  point  for  doublet 
lattice 

pressure  (acting  at  1/4  chord  for  doublet  lattice) 
aerodynamic  Influence  matrix. 

(Ref.  6:17.5-1) 


9 


where 


Ajj(k,m)  - 


NASTRAN  takes  all  codes  for  computing  A^^  directly  from  Gleslng,  Kalman, 
and  Rodden  (Ref.  A).  The  dovnwash  at  the  3/4  chord  point  Is  related  to 
the  deflection  and  rotation  of  the  box  control  point  at  the  center.  After 
solving  for  the  unknown  pressure  distribution  of  the  wing  from  Eq.  (2),  a 
spline  fitting  technique  Is  used  to  relate  the  set  of  point  aerodynamic 
loads  to  the  structural  grid. 

Surface  Spline 

The  structural  grid  points  usually  do  not  coincide  v/lth  the  grid 
points  that  define  the  aerodynamic  elements.  Equations  of  restraint 
between  the  two  sets  of  grid  points  are  supplied  by  the  theory  of  surface 
splines  (Ref.  6).  A surface  spline  Is  a solution  for  an  Infinite  uniform 
plate,  or  a function  h(x,y)  for  all  points  (x,y),  when  a discrete  set  of 
points,  h^  •>  h(x^,y^).  Is  known.  In  NASTRAN,  the  structural  degrees  of 
freedom  are  taken  to  be  the  Independent  degrees  of  freedom,  and  the  aero- 
dynamic degrees  of  freedom  are  dependent.  From  the  theory  of  surface 
splines,  then,  a matrix  Is  obtained  which  relates  the  dependent  degrees 
of  freedom  to  the  Independent  degrees  of  freedom.  This  feature  of 
NASTRAN  allows  the  choice  of  structural  and  aerodynamic  elements  to  be 
made  Independently  of  each  other. 

Flutter  Equation 

The  basic  differential  equation  for  an  oscillating  lifting  surface 
Is 


[M]  {U}  - [A]  {u}  + [K]  {u}  - 0 


(3) 


where 
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[M]  > generalized  mass  matrix 

[A]  « generalized  aerodynamic  force  matrix 

[K]  » generalized  stiffness  matrix 

{u}  ■ generalized  modal  vector. 

(Ref.  15:92) 

The  generalized  aerodvnamlc  force  matrix  Is  the  product  of  the  djmamlc 
2 

pressure,  1/2  pU  , and  the  aerodynamic  matrix,  [Q],  where  an  element  of 
[Q]  Is  defined  as 


(4) 


where  Q (x,y)  Is  the  deflection  distribution  of  the  r^^  mode,  C , ,(x,y) 
r pCs) 

Is  the  pressure  coefficient  distribution  In  the  s^^  mode,  and  the  Inte- 
gration Is  performed  over  the  lifting  surface  area  S.  Then,  Equation  (3) 
may  be  written  as 

,2 


[M]  (u)  - [Q]  {u}  + [K]  {u}  - 0 


(5)  . 


To  Include  structural  damping  In  the  system,  the  generalized  stiffness 
matrix  becomes 


[K]  - [I  + 1 G]  [K] 


(6) 


where  [G]  Is  a diagonal  matrix  of  modal  damping  and  [I]  Is  an  Identity 
matrix.  The  stiffness  In  the  mode  Is  then 


(1  + 


jJ 


(7) 


where  g^^^^ 
(5)  becomes 


Is  the  structural  damping  In  the  j 


th 


mode. 


Then  Equation 
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[M]  {u}  - ^ [Q]  {u}  + [ic]  {u}  = 0 


Assume  simple  harmonic  motion.  Then  the  generalized  coordinates  are 
written  as  {u}  = {u}  where  u is  no  longer  time  dependent  and 

may  contain  damping.  Equation  (8)  is  then  rewritten  as 


[^2  _ P|_  [Q]  + [K]]  {u}  = 0 


This  is  the  basic  flutter  equation. 


K-Method  of  Flutter  Solution 

In  the  k-method  of  flutter  solution,  fi  is  considered  to  be  undamped 
(£5  “ lu) , and  the  aerodynamic  force  matrix  is  a function  of  reduced  fre- 
quency k and  Mach  number  m.  Then,  the  basic  equation  for  modal  flutter 


analysis  is 


[Q(k,m)]  + [K]]  {u}  = 0 


(Ref.  15:93) 


The  reduced  frequency  is  defined  as 


, COJ 

“ 2U 


where 

c “ a reference  length 
u ••  the  frequency  of  oscillation 
U “ the  free  stream  velocity. 

In  order  to  bring  the  system  into  neutral  stability,  an  artificial  struc- 
tural damping  term,  g,  is  introduced.  Applying  this  artificial  damping 
to  the  diagonal  term  of  the  stiffness  matrix  of  Equation  (10)  yields 


( 


Kjj  - ^ (1  + Igg(j)  + Ig)  Kjj  (12) 

For  small  values  of  damping  (g*g  (j)<<l)*  the  diagonal  term  is 

s 

approximated  by 

Kj.  - (^)  (1  + Igg(j))  K (13) 


or 


^ 2 ^ 
GJ 


The  aerodynamic  terms  are  then  converted  to  aerodynamic  mass : 


- [M]  + I (-^)Mq  (k.m)] 


1+ig 

I- 


+ [K] 


{u}  = 0 


(lA) 


i 

(15) 


Equation  (15)  is  a complex  eigenvalue  problem  that  is  solved  for  a series 

of  values  for  parameters  k,  m,  and  p.  The  complex  eigenvalues  are 
2 

u /(1+lg),  which  yield  real  values  of  circular  frequency  o)  and  structural 
damping,  g.  Velocity  U is  obtained  from  U = c(jj/2k.  A positive  value  of 
damping  requires  structural  damping  be  added  to  bring  the  system  into 
neutral  stability,  or,  the  system  may  be  considered  unstable.  A negative 
value  has  the  opposite  interpretation.  Flutter  occurs  when  the  values  of 
k,  m,  and  p force  the  artificial  structural  damping,  g,  to  be  zero. 

NASTRAN  uses  a variation  of  Eq.  (15)  that  allows  solution  of  static 
divergence  problems.  The  equation  is 

[[(^)^[M]  + (f)  “ 0 (16) 

(Ref.  7:17.6-1) 
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L 


■'jj  ■ 7 ■'jj 


(12) 


.th 


For  small  values  of  damping  (g’g^Cj )<<1) , the  J diagonal  term  Is 


approximated  by 


(13) 


or 


ir  " M TT 

ii  ^ ^2  ^ ^jj 


(14) 


The  aerodynamic  terms  are  then  converted  to  aerodynamic  mass : 


- [M]  + f (^)MQ  (k.m)] 


L 


r 2 

u 


1+lg 


L ^ 


+ [ic] 


] (u)  - 0 


(15) 


Equation  (15)  Is  a complex  eigenvalue  problem  that  Is  solved  for  a series 


of  values  for  parameters  k,  m,  and  p.  The  complex  eigenvalues  are 
2 


0)  /(1+lg),  which  yield  real  values  of  circular  frequency  ui  and  structural 
damping,  g.  Velocity  U Is  obtained  from  U “ ca)/2k.  A positive  value  of 
damping  requires  structural  damping  be  added  to  bring  the  system  Into 
neutral  stability,  or,  the  system  may  be  considered  unstable.  A negative 
value  has  the  opposite  Interpretation.  Flutter  occurs  when  the  values  of 
k,  m,  and  p force  the  artificial  structural  damping,  g,  to  be  zero. 

NASTRAN  uses  a variation  of  Eq.  (15)  that  allows  solution  of  static 
divergence  problems.  The  equation  Is 


[[(f)^[M]  + (f)  tQ(k.m)i][^]+  [ici]  {G}  - 


(16) 


(Ref.  7:17.6-1) 
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It  is  then  possible  to  let  k ••  0,  as  division  by  k Is  not  required.  At 
k ■ 0,  w ■ 0,  which  indicates  the  static  condition  of  divergence  exists 
if  damping  becomes  positive. 


N 1 

’ 4 

i III.  Analysis 

I 

I 


Structure  Description 

The  wing  box  top  view  Is  shown  In  Fig.  3 with  Its  node  points 

numbered  as  shown.  The  corresponding  lower  node  point  numbers  are 

found  by  adding  one  to  the  top  node.  The  element  gage  sizes  are  shown 

In  Fig.  A,  where  the  bottom  plate  elements  have  the  same  dimensions  as 

the  top  plate  elements.  The  rod  elements  that  connect  the  upper  and 

2 

lower  nodes  all  have  a cross  sectional  area  of  0.10  In  . The  wing 
dimensions  and  element  sizes  are  the  same  as  the  strength  optimized 
Intermediate  complexity  wing  from  AFFDL-TR-75-137  (Ref.  15:167).  Fig.  5 
gives  the  wing  box  dimensions. 

Elements.  The  wing  box  consists  of  100  finite  elements.  The  upper 
and  lower  surfaces  of  the  wing  box  are  covered  by  trapezoidal  plate 

i 

elements  ~ 20  on  the  upper  surface  and  20  on  the  lower  surface.  A total  j 

of  35  shear  panels  represent  the  three  spars  and  five  ribs.  To  Inhibit  ! 

motion  between  the  upper  and  lower  surfaces,  25  rods  connect  the  upper  ] 

and  lower  surface  node  points.  | 

( 

Degrees  of  Freedom  at  Nodes.  There  are  a total  of  60  nodes  at 
which  the  100  elements  connect.  The  ten  nodes  In  the  root  chord  plane 
are  rigidly  constrained  In  order  to  cantilever  the  wing.  Each  of  the 
remaining  50  nodes  are  allowed  one  degree  of  freedom  normal  to  the  plane 
of  the  wing. 

In  addition,  the  nodes  along  the  leading  and  trailing  edges  are 
allowed  a rotational  degree  of  freedom  about  axes  running  along  the 
leading  and  trailing  edge.  The  tip  nodes  are  allowed  to  have  two  rota- 
tional degrees  of  freedom,  as  shown  In  Fig.  6.  These  additional  degrees 
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of  f/eedom  are  necessary  in  order  to  simulate  the  rotational  Inertial 
, effects  of  the  overhanging  panels  connected  to  the  wing  box. 

Wing  Mass.  The  fully  stressed  wing  structure  weight  was  Initially 
61.1  lbs.  Then,  additional  mass  was  added  to  account  for  structural 
connection  materials,  such  as  flanges  and  rivets,  and  nonstructural  mass  i 

such  as  pumps,  tubing,  etc.  In  addition,  weight  for  the  overhanging 
structure  was  Incorporated  at  the  leading  edge,  trailing  edge,  and  tip 
nodes,  with  representative  moments  of  Inertia  about  these  nodes  to 

simulate  rotatlotal  Inertial  effects  of  the  panels.  The  lumped  masses  j 

and  moments  at  each  node  are  found  In  Appendix  A.  These  additional 
masses  raised  the  total  weight  of  the  wing  to  374.6  lbs. 

Aerodynamic  Description 

The  wing  flutter  analysis  was  performed  using  the  doublet-lattice 


aerodynamics  for  Mach  0.9  at  sea  level.  The  planform  was  represented 
by  32  trapezoidal  aerodynamic  panels,  four  chordwlse  and  eight  spanwlse 
at  equal  division  points,  as  seen  In  Fig.  7.  These  divisions  were  chosen 
In  order  to  make  the  aspect  ratio  of  the  boxes  roughly  one  or  less.  Also, 
the  number  of  boxes  per  reference  chord  should  be  greater  than  eight 
times  the  reduced  frequency  (Ref.  8:1.11-1).  The  aerodynamic  model 
does  not  Include  the  effect  external  stores  would  have  upon  the  lifting 
surface. 

A surface  spline  was  used  to  Interpolate  between  all  the  aerodynamic 
grid  points,  located  at  the  center  of  each  aerodynamic  element,  and  the 
selected  structural  grid  points.  The  odd  numbered  structural  nodes  were 
selected  for  spline  Interpolation  with  all  the  aerodynamic  nodes.  Only 
the  normal  displacement  degree  of  freedom  for  each  structural  node  Is 
required  for  the  surface  spline. 


i 


20 


I 


Elastic  Axis  Determination 

The  location  of  the  elastic  axis  of  the  vlng  was  essential  in  order 
to  correctly  Interpret  the  flutter  results.  The  elastic  axis  Is  defined 
as  the  axis  about  which  rotation  occurs  when  a wing  Is  loaded  In  pure 
torsion  (Ref.  9:448).  It  may  alternately  be  defined  as  the  locus  of  shear 
centers  of  the  cross  sections  of  a cantilever  wing  (Ref.  3:17).  The  shear 
center  of  a wing  cross  section  Is  the  point  at  which  a shear  force  acts  to 
produce  a wing  deflection  without  rotation.  If  the  wing  is  an  elastic 
structure,  the  shear  center  of  a cross  section  coincides  with  the  elastic 
axis. 

To  locate  the  elastic  axis  on  the  wing  model,  a static  analysis  was 
performed  using  NASTRAN.  A static  shear  force  was  applied  to  each  lower 
node  point  along  each  rib.  The  elastic  axis  was  located  by  interpolation 
between  the  two  nodes  that  produced  the  least  rotation  of  the  cross  sec- 
tion. The  resulting  locus  of  shear  centers  Is  shown  In  Fig.  8. 

Free  Vibration 

A free  vibration  analysis  was  performed  on  the  basic  wing  model  to 
determine  the  fundamental  mode  shapes  and  frequencies.  Only  vertical 
displacement  of  nodes  11  through  60  were  used  in  the  analysis  with  nodes 
1 through  10  being  totally  restrained.  The  first  three  fundamental  modes 
obtained  are  seen  In  Figs.  9 through  11,  In  which  the  deformed  shape  Is 
superimposed  over  the  undeformed  shape.  Mode  1 Is  seen  to  be  first  bend- 
ing with  a frequency  of  49.6  Hz.  Mode  2 Is  predominantly  first  torsion 
coupled  with  some  bending.  This  frequency  is  95.4  Hz.  Mode  3 Is  primarily 
second  bending  coupled  with  considerable  torsion.  This  frequency  Is  104.1 
Hz.  These  three  modes,  used  for  the  modal  flutter  analysis,  are  typical 
for  a simple  cantilevered  swept  wing. 
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Location  of  Elastic  Axis 


I Flutter  Analysis 

All  flutter  analysis  is  performed  for  conditions  at  sea  level  and 
Mach  0.9.  The  first  flutter  test  was  for  the  basic  wing  model  with  no 
external  stores  added.  A reduced  frequency  range  from  k=0.75  to  k“0.0 
was  used  to  check  for  flutter.  At  k=0.0,  flutter  frequency  goes  to  zero, 
which  indicates  the  static  condition  of  divergence  occurs  if  damping 
becomes  positive  at  that  point.  As  mentioned  in  Chapter  II,  NASTRAN 
uses  a modified  flutter  equation,  Eq.  (16),  which  allows  k to  go  to  zero. 

The  next  flutter  test  determined  the  effect  of  concentrated  masses 
placed  along  ribs  and  spars  at  nodes  of  the  wing.  This  procedure  is 
commonly  called  mass  balancing.  Two  mass  values  were  used  for  the  test. 

A 3.10  slug  mass  (100  lb),  or  approximately  27%  of  the  total  wing  mass, 
and  a 6.21  slug  mass  (200  lb),  or  approximately  53%  of  the  total  wing 
mass,  were  selected.  For  each  node  that  received  a mass,  a flutter  speed 
and  flutter  frequency  were  determined. 

To  simulate  store  pitch  inertia,  two  equal  masses  were  used.  In 
order  to  compare  results  with  the  single  concentrated  mass  results,  the 
centers  of  gravity  and  mass  value  for  both  cases  were  the  same.  As  an 
example,  to  simulate  a 100  lb  store  with  pitch  Inertia,  two  50  lb  weights 
were  placed  along  the  rib  on  the  two  nodes  adjacent  to  where  the  single 
concentra'.ed  ]00  lb  weight  had  been.  The  same  procedure  was  followed  for 
the  2C0  lb  store  to  simulate  store  pitch  inertia.  Fo  each  mass  combina- 
tion, a flutter  speed  and  flutter  frequency  were  determined. 

The  normal  procedure  to  determine  flutter  characteristics  was  to 
first  use  a wide  range  of  k values  to  find  two  k values  where  the  value 
of  artificial  damping  changed  from  negative  to  positive.  Then,  a range 
of  k values  between  the  two  k values  was  selected  to  closely  bracket  the 

s 
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flutter  velocity  and  flutter  frequency  values  where  damping  goes  to  zero. 

A linear  variation  of  velocity  and  frequency  was  then  assumed  for  the  In- 
terpolation to  the  values  of  velocity  and  frequency  that  occur  at  zero 
damping.  This  assumption  Is  accurate  If  the  flutter  frequency  and  velocity 
are  bracketed  closely. 

For  problems  that  require  considerable  amounts  of  computer  time,  a 
saving  In  time  could  be  realized  when  refining  the  location  of  the  flutter 
speed  by  using  the  restart  procedure  Incorporated  In  NASTRAN.  By  using 
the  restart  procedure,  recomputation  of  the  natural  modes  and  frequencies 
of  the  structure  would  be  avoided.  Likewise,  only  newly  requested  aero- 
dynamic matrices  would  be  computed,  which  Is  a sizeable  saving  In  computer 
time. 
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IV.  Results 


The  two  flutter  plots  generated  for  each  case  are  represented  by 
example  Figs.  12-A  and  12-B.  The  flutter  plots  for  each  case  in  which 
flutter  occurred  are  found  in  Appendix  B.  Figure  12-A  is  the  damping 
versus  velocity  plot,  from  which  the  flutter  speed  and  flutter  mode  may 
be  determined.  Figure  12-B  is  the  frequency  versus  velocity  plot,  from 
which  the  flutter  frequency  is  determined  by  knowing  the  flutter  velocity 
and  flutter  mode. 

The  values  for  flutter  speed  and  frequency  obtained  by  linear  inter- 
polation for  the  various  cases  are  listed  in  Tables  I and  II.  Table  I 
contains  the  values  obtained  for  the  mass  balancing  experiment  where  single 
concentrated  100  lb  and  200  lb  weights  were  placed  at  various  nodes. 

Table  II  lists  the  flutter  speeds  and  frequencies  obtained  for  the  store 
inertia  simulation,  plus  the  computed  store  pitch  inertia  about  the 
store  center  of  gravity.  The  investigative  Mach  number  was  0.9,  or  577 
knots  for  conditions  at  sea  level,  for  all  cases.  The  flutter  speeds 
obtained  were  all  much  greater  than  577  knots,  which  means  that  this 
wing  model  will  not  flutter  at  sea  level  at  0.9  Mach.  This  is  a desirable 
characteristic. 

For  the  basic  wing  with  no  external  stores  added,  no  flutter  was 
detected  up  to  divergence  speed.  The  wing  diverged  in  the  torsion  mode. 
This  result  indicates  a wing  that  has  sufficient  torsional  stiffness  to 
prevent  flutter  up  to  wing  divergence  when  it  carries  no  external  stores. 

The  free  vibration  frequencies  for  each  loading  case  are  given  in 
Table  IV  of  Appendix  C.  Also  in  Appendix  C,  plots  comparing  the  free 
vibration  frequencies  of  store  versus  concentrated  mass  at  the  3/4  node 
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Fig.  12-B.  Flutter  Frequency  Solution 


Table  I 

Results  for  Mass  Balancing 


100 

lb 

Concentrated  at 

a Node 

Node 

(knots) 

(Uj  (Hz) 

28 

none 

__  _ 

30 

none 

— 

36 

none 

— 

38 

2357 

45.6 

40 

2096 

49.4 

44 

none 

46 

2226 

44.1 

48 

2169 

41.6 

50 

1961 

48.0 

54 

none 

... 

56 

2133 

39.8 

58 

2136 

36.7 

60 

1937 

42.2 

200 

lb 

Concentrated  at 

a Node 

28 

2262 

42.5 

38 

2387 

38.1 

48 

2243 

33.5 

58 

2182 

28.8 
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Table  II 

Results  for  Store  Inertia 

Simulation 

100  lb  Store  (Two  x 

50  lb) 

Nodes 

CG  Node 

*Icg(lb-ft-sec^) 

(knots) 

0)^  (Hz) 

26, 

30 

28 

342 

none 

34, 

38 

36 

277 

none 

— 

36, 

40 

38 

277 

2414 

49.5 

44, 

48 

46 

219 

none 

— 

46, 

50 

48 

219 

2153 

48.9 

54, 

58 

56 

175 

2618 

36.7 

56, 

60 

58 

175 

2139 

50.7 

200  lb  Store  (Two  x 

100  lb) 

26, 

30 

28 

685 

2435 

49.8 

36. 

40 

38 

555 

2263 

48.4 

46, 

50 

48 

438 

2130 

44.5 

56, 

60 

58 

350 

2270 

40.4 

*Note:  leg  Is  store  pitch  Inertia  about  the  store  center  of 
gravity 


points  are  given.  The  mode  shapes  for  each  loading  case  were  not  obtained. 


Mass  Balancing  Effects 

Masses  were  added  to  the  basic  wing  using  one  mass  value  at  a node 
for  each  flutter  test.  The  mass  was  shifted  in  both  the  spanwlse  and 
chordwlse  directions.  The  first  observation  was  that  no  flutter  could  be 
produced  by  placing  the  mass  between  the  leading  edge  and  the  elastic 
axis.  This  is  a normal  result  as  flutter  Is  damped  In  most  cases  where 
the  center  of  gravity  of  a wing  is  shifted  forward  of  the  elastic  axis. 
When  the  mass  was  shifted  aft  of  the  elastic  axis,  flutter  was  produced 
except  for  those  cases  of  the  100  lb  weight  placed  near  the  wing  root. 

As  can  be  seen  In  Table  I,  as  the  mass  was  shifted  chordwlse  from 
the  nodes  at  midchord  along  a rib  to  the  nodes  at  the  trailing  edge,  the 
flutter  speed  generally  decreased.  The  lowest  flutter  speeds  occurred 
when  the  mass  was  placed  on  the  trailing  edge  nodes  of  each  rib.  The 
most  critical  flutter  speed  was  produced  by  placing  the  mass  at  node  60 
located  at  the  tip.  No  definite  trend  for  flutter  frequency  could  be 
detected- for  the  chordwlse  movement  of  mass.  As  the  mass  was  moved 
spanwlse  toward  the  tip  of  the  wing,  the  flutter  speed  and  frequency 
both  decreased.  When  moving  the  mass  spanwlse,  the  critical  flutter 
speed  occurred  when  the  mass  was  located  at  the  tip  of  the  wing. 

The  flutter  mode  was  the  fundamental  torsion  mode  In  all  cases. 
Evident  on  most  of  the  frequency  versus  velocity  plots  in  Appendix  B Is 
the  Interaction  of  the  first  two  modes  when  flutter  occurs.  The  Inter- 
action Is  characterized  by  the  frequencies  of  the  two  modes  converging, 
or  nearly  so,  to  a single  frequency  at  flutter.  This  Is  the  classic 
coalescence  of  frequencies  that  Is  a characteristic  of  modal  flutter. 
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Fig*  13.  Comparison  of  Flutter  Speeds  of  100  lb  Store  Versus 
Single  Concentrated  Mass  at  3/4  Chord  Nodes 


Store  Inertia  Effects 


The  results  of  the  flutter  test  Involving  two  50  lb  weights,  simulat- 
ing a 100  lb  store,  are  compared  with  results  for  the  single  concentrated 
100  lb  weight  at  the  3/4  chord  nodes  in  Figs.  13  and  14.  Looking  at 
Fig.  13  first,  the  largest  difference  in  flutter  speed  for  the  two  cases 
was  at  node  38.  Here,  the  store  pitch  inertia  raised  the  flutter  speed 
over  the  single  concentrated  mass  flutter  speed  by  fifty-seven  knots. 

There  is  a definite  trend  in  Fig.  14,  as  the  flutter  frequency  of  the 
store  simulation  is  consistently  higher  than  the  single  mass  case.  The 
difference  increases  from  3.9  Hz  at  node  38  to  14.0  Hz  at  node  58. 

Now  consider  Figs.  15  and  16,  and  compare  flutter  results  of  the  200 
lb  store  (two-mass  combination)  with  the  single  200  lb  mass  along  the  3/4 
chord  line.  As  can  be  seen  in  Fig.  15,  the  two-mass  combination  tends  to 
damp  flutter  at  nodes  38  and  48,  and  conversely,  to  promote  flutter  at 
nodes  28  and  58,  relative  to  the  single  mass.  Thus,  there  is  no  definite 
trend  with  flutter  speeds  as  affected  by  store  pitch  inertia.  But  look- 
ing at  flutter  frequencies  in  Fig.  16,  there  is  again  a definite  trend 
as  the  flutter  frequency  of  the  simulated  store  is  consistently  higher 
than  the  single  mass  flutter  frequency.  The  difference  varied  from  7.3 
Hz  at  ncde  28  to  11.6  Hz  at  node  58. 
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V,  Conclusions  and  Recommendations  for  Further  Study 


Conclusions 

The  conclusion  that  can  be  drawn  from  the  mass  balancing  analysis 
is  that,  as  the  center  of  gravity  of  the  wing  moves  aft  of  the  elastic 
axis,  the  flutter  speed  of  the  wing  is  reduced.  Likewise,  as  mass  is 
shifted  spanwlse  toward  the  tip,  the  flutter  speed  and  frequency  are 
both  reduced.  It  should  not  be  inferred  that  all  swept  back  wings  will 
react  in  this  manner,  as  the  mode  shape  interactions  differ  from  wing  to 
wing  due  to  mass  distribution  and  stiffness  variations.  For  example,  in 
the  experimental  flutter  tests  of  Runyan  and  Sewall  (Ref.  11:6),  the 
flutter  modes  were  first  and  second  bending,  versus  first  torsion  for 
the  wing  used  in  this  study.  Nevertheless,  mass  balancing  produced 
definite  flutter  speed  trends  for  this  particular  wing. 

In  the  external  store  simulation,  there  was  a definite  increase  in 
the  flutter  frequency  due  to  store  pitch  inertia  when  compared  to  the 
flutter  frequency  of  a concentrated  mass  at  the  store  center  of  gravity. 

No  definite  trend  was  established  for  flutter  speed,  however.  The  upward 
shift  in  flutter  frequency  indicated  that  the  torsional  mode  had  a 
stronger  participation  in  the  simulation  than  with  the  single  concentrated 
mass. 

As  for  the  method  of  analysis,  NASTRAN  is  a versatile  structural 
analysis  computer  program  and  is  relatively  easy  to  use.  However,  the 
user  must  have  a good  understanding  of  finite  elements,  and  now  must  be 
versed  in  aerodynamic  theory  as  well  to  successfully  use  the  flutter 
format.  To  add  versatility  to  the  flutter  package,  the  p-k  method 
needs  to  be  added,  as  it  gives  better  subcritlcal  frequency  and  damping 
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trends  than  the  k method.  Also,  to  supplement  the  subsonic  doublet-lattice 


aerodynamics,  supersonic  capability  should  be  added  to  NASTRAN  by  including 
Mach  Box  or  Kernel  Function  aerodynamic  theory  in  the  program. 

Recommendations  for  Further  Study 

For  further  stores  flutter  clearance  studies,  a wing  model  of  an 
operational  aircraft,  such  as  the  F-16,  should  be  considered  so  that  a 
comparison  with  experimental  results  is  possible.  Another  possibility 
is  the  modeling  of  store-wing  aerodynamic  interaction  to  see  what  effect 
there  is  on  flutter,  versus  a model  with  no  Interaction,  as  in  this 
study.  Finally,  a more  accurate  mass  model  of  a store  could  be  used. 

For  example,  a store  that  is  flexibly  connected  to  the  wing  structure 
would  be  a more  realistic  representation  of  store-wing  dynamic  inter- 
iut  ion. 
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Appendix  A 


Table  III 
Mass  and  Mass  Moments  of 

Inertia  at  Nodes 

Nodes 

Mass  at  Each  Node 
(slugs) 

Mass  Moment  of  Inertia 
(slug-ft^) 

11, 

12 

0.6525 

0.0324 

13, 

14 

0.1798 

15, 

16 

0.2919 

17, 

18 

0.2869 

19, 

20 

0.4438 

0.0259 

21, 

22 

0.4177 

0.0125 

23, 

24 

0.1873 

25, 

26 

0.2633 

27, 

28 

0.2345 

29, 

30 

0.3444 

0.01725 

31, 

32 

0.4000 

0.0097 

33, 

34 

0.1627 

35, 

36 

0.2093 

37, 

38 

0.1696 

39, 

40. 

0.3028 

0.0140 

41, 

42 

0.2484 

0.0054 

43, 

44 

0.0928 

45, 

46 

0.1106 

47, 

48 

0.0792 

49, 

50 

0.2174 

0.0151 

51, 

52 

0.1437 

0.0058 

*53, 

54 

0.0838 

(1)  0.0058;  (2)  0.0049 

55, 

56 

0.0937 

*57, 

58 

0.0758 

(1)  0.0051;  (2)  0.0051 

59, 

60 

0.1242 

0.0086 

*Note:  All  moments  are  about  a local  axis  parallel  to  edge  of 
structure  except  at  wing  tip  where  second  value  of 
Inertia  is  about  the  normal  axis  depicted  in  Fig.  6. 


Appendix  B 


Flutter  Speed  and  Flutter  Frequency  Graphs 
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Fig.  18-B.  Flutter  Frequency  Solution  with  100  lb  at  Node  46 
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Fig.  19-A.  Flutter  Speed  Solution  with  100  lb  at  Node  48 
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19-B.  Flutter  Frequency  Solution  with  100  lb  at  Node  48 


Fig.  20-B.  Flutter  Frequency  Solution  with  100  lb  at  Node  50 
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Flutter  Frequency  Solution  with  100  lb  at  Node  56 
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Fig.  22-B.  Flutter  Frequency  Solution  with  100  lb  at  Node  58 
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Flutter  Frequency  Solution  with  50  lb  at  Node  36 
and  50  lb  at  Node  40 


Fig.  25-A.  Flutter  Speed  Solution  with  50  lb  at  Node  46 
and  50  lb  at  Node  50 
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Fig.  25-B.  Flutter  Frequency  Solution  with  50  lb  at  Node  46 
and  50  lb  at  Node  50 
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Fig.  26-B.  Flutter  Frequency  Solution  with  50  lb  at  Node  64 
and  50  lb  at  Node  58 
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Fig.  32-B.  Flutter  Frequency  Solution  with  100  lb  at  Node  26 
and  100  lb  at  Node  30 
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Pig.  35-B.  Flutter  Frequency  Solution  with  100  lb  at  Node  56 
and  100  lb  at  Node  60 
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Fig.  36.  Coaparlson  of  Mode  1 Free  Vibration  Frequencies  of  100  lb 
Store  Versua  Single  Concentrated  Mass  at  3/4  Chord  Nodes 
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Fig.  37.  Comparison  of  Mode  2 Free  Vibration  Frequencies  of  100  lb 
Store  Versus  Single  Concentrated  Mass  at  3/4  Chord  Nodes 
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Fig.  38.  Comparison  of  Mode  1 Free  Vibration  Frequencies  of  200  lb 
Store  Versus  Single  Concentrated  Mass  at  3/4  Chord  Nodes 
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